Energetics and Possible formation and decay mechanisms of Vortices in Helium 

Nanodroplets 
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The energy and angular momentum of both straight and curved vortex states of a helium nano- 
droplet are examined as a function of droplet size. For droplets in the size range of many experiments, 
it is found that during the pickup of heavy solutes, a significant fraction of events deposit sufficient 
energy and angular momentum to form a straight vortex line. Curved vortex lines exist down to 
nearly zero angular momentum and energy, and thus could in principle form in almost any collision. 
Further, the coalescence of smaller droplets during the cooling by expansion could also deposit suffi- 
cient angular momentum to form vortex lines. Despite their high energy, most vortices are predicted 
to be stable at the final temperature (0.38 K) of helium nanodroplets due to lack of decay channels 
that conserve both energy and angular momentum. 



Vortices are an almost unavoidable presence in bulk 
superfluid helium. There is a rich history of studies of 
their properties and interactions X It is natural to con- 
sider their possible presence in the finite superfluid found 
in 4 He nanodroplets^ The study of helium nanodroplets 
has been quite active in recent years, starting with the 
spectroscopy of embedded molecules as a probe for prop- 
erties of nanodroplets, which was introduced by Scoles 
and coworkers^ However, despite a now large body of 
work, no unambiguous signature of the presence of vor- 
tices has yet been reported. This is perhaps surprising, 
as a vortex in a nanodroplet is expected to bind a molec- 
ular impuritji^ and likely introduce a highly anisotropic 
interaction potential, though no explicit calculation of 
the magnitude of such an anisotropy has been reported 
to date. The vortex-induced anisotropy in molecular ori- 
entation, if large compared to the rotational constant of 
the molecule in liquid helium, will quench the molecu- 
lar rotation and collapse the rotational structure that is 
one of the hallmarks of molecular spectroscopy in helium 
nanodropletsi&i 

Several calculations have been published for the en- 
ergy of a droplet with a straight vortex, with and with- 
out an atom or cylindrically symmetric molecule aligned 
with the vortexi 4 *^*^ As in bulk superfluid, these calcu- 
lations have found the energy of droplets with a straight 
vortex to be significantly higher than that of vortex-free 
droplets4£ This result, combined with the failure to ob- 
serve vortex lines to date, have led some to propose that 
vortices are unstable in helium nanodroplets and per- 
haps are rapidly expelled. In this paper, we revisit the 
energetics of vortex lines, considering both the straight 
vortex line down the center of the droplet that has pre- 
viously been considered and curved vortex lines that will 
rotate around the droplet due to their own flow field. 
Furthermore, we consider the effect of angular momen- 
tum conservation on the possible formation and decay 
of vortices in droplets. It turns out that the linear and 
curved vortices contain considerably less excitation en- 
ergy per unit angular momentum than the final states 
accessible by decay through the mechanisms considered. 
This suggests that despite their higher energy than the 



ground state, droplets with such vortices are the lowest 
possible states with high angular momentum, and thus 
should in fact be stable to decay. 



I. HOLLOW CORE MODEL OF VORTEX LINES 

There have been several microscopic calculations of 
the properties of straight vortex lines in helium nan- 
odroplets, both pure and doped with atomic or molec- 
ular solutesA^i&iLiSiii In this work, we exploit the phe- 
nomenological description known as the hollow core 
model. The numerous microscopic treatments of vortex 
lines in two and three dimensions have largely confirmed 
its qualitative applicability^SiiiiliiSiiSiii and it is widely 
used to describe the properties of vortex lines and rings 
in bulk liquid hcliumi 1 ^ 

The vortex is surrounded by a circulating flow, char- 
acterized by an irrotational velocity field (V x v = 0). 
In the simple case of a straight vortex in bulk helium, 
the magnitude v of the flow velocity at any given point 
is inversely proportional to the distance of that point 
from the center of the vortex, v = h/mr (m is the mass 
of a 4 He atom). This leads to a kinetic energy den- 
sity whose volume integral diverges as the vortex is ap- 
proached (r — > 0). In the hollow core model, the helium 
number density p is taken to be zero inside a cylinder of 
radius a and equal to the bulk value pb = 0.0218 A~ 3 out- 
side of this cylinder. Such a discontinuous change in den- 
sity is unphysical, but microcanonical calculations have 
confirmed a nearly hollow core, though with a smooth 
transition of the density to the bulk valueiii 1 ^ It is noted 
that our model treats the helium density on the outer 
boundary of the droplet as abruptly going to zero, while 
it is known that in fact the surface of liquid helium is 
diffuse, with a thickness of w 6 — 8 A. 19 A hollow core 
radius of a = 1.00 A was found experimentally*^ to best 
reproduce the measured energy and velocity of vortex 
rings in bulk liquid helium; matching our expression of 
the velocity of vortex rings to the expression in Ref. Il8l 
requires us to use a = a/y/e — 0.607 A (see Appendix lA")l . 
Any normal fluid component is neglected as helium nan- 
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odroplets in this size range have no thermally excited 
phonon excitations S 

First we consider a straight vortex through the center 
of a helium nanodroplet. Because of the loss of spherical 
symmetry, we will consider that droplet to have the shape 
of an ellipsoid, with axial radius B and equatorial radius 
A; the density is assumed to be uniform, equal to the 
bulk value pb, and dropping to zero at the ellipsoidal 
surface. In this geometry, the flow field is the same as 
for the bulk straight vortex described above. The kinetic 
energy, E sv , of this straight vortex flow field is given by 



E„ 




VA 2 



A 



for A > a. (1) 



For a spherical droplet of N atoms, B = A = R = 
^/r 2 A 2 / 3 + a 2 , with r = (47rp b /3) _1/3 = 2.22 A. 

An interesting question is the extent to which the large 
angular momentum of the vortex will distort the other- 
wise spherical droplet. Making the droplet oblate will 
reduce the length of the vortex line and thus lower its en- 
ergy. However, this will also increase the surface energy 
of the droplet. The latter is the product of the surface 
tension of bulk liquid helium, a = 0.272 K A~ 2 , — and the 

surface area of an ellipsoid, S = 2irA 2 + 



■ln(i±£ 

1 1 — e 



with the eccentricity e = \J\ — (B/A) 2 . The droplet 
distortion is found by minimizing the sum of the vor- 
tex kinetic energy (given by Eq. ijTjl) and the surface 
energy of the droplet, at constant droplet volume V = 
^-^(A 2 — a 2 ) 3 / 2 . Table|Jgives the resulting eccentricity, 
e, and the associated stabilization energy (reduction in 
energy from a spherical droplet of the same volume with 
vortex). It is seen that the distortion from spherical sym- 
metry is small, in the sense that B/A = \f\ — e 2 is close 
to 1, despite the high angular momentum of the vortex, 
and this distortion will be neglected in the rest of this 
paper. 

Dalfovo et alA used finite range Density Functional 
Theory (DFT) to calculate the energy and core shape of 
a straight line vortex for a range of droplet sizes (N = 
50—1000). While their core has a smooth density profile, 
they report a core radius of the order of 1 — 2 A. Figure^] 
shows a plot of the vortex energy calculated by the hollow 
core model as a function of the number of helium atoms 
in the droplet, along with the same quantity estimated by 
the Density Functional 4 method. The hollow core model 
appears to slightly underestimate the vortex energy with 
respect to DFT; agreement with the DFT calculation can 
be made almost quantitative if a value of a — 0.56 A is 
used. The core radius used in the hollow core model is 
to be interpreted as an effective radius that reproduces 
the energetics and velocity of experimental vortices, and 
not more than an estimate of the real core radius. 

The angular momentum associated with a straight vor- 
tex is NH, i.e., one unit per helium atom. Both methods 




FIG. 1: Energy of the straight vortex line solution as a func- 
tion of the number of He atoms in the droplet, calculated 
both by the Hollow Core model (solid line, Eq. and by 
finite range Density Functional Theory (dashed line, Ref. 0). 
Energies are displayed as a fraction of the minimum energy 
required to deposit NTi units of angular momentum into rip- 
plon excitation modes. 



agree that the vortex energy per unit angular momentum 
is lower than for any other excitation mode of a pure 
helium droplet: 20 Figure ^ shows the vortex excitation 
energies relative to the energy required to produce N /2 
quanta of L — 2 ripplons, which is the lowest energy state 
based upon quasiparticle excitations that has the same 
total angular momentum as the straight vortex. Table |I] 
shows a comparison of the hollow core model excitation 
energy of a straight vortex and this lowest ripplon, and 
their energies per unit angular momentum for droplets 
of the size range that span most helium nanodroplet ex- 
periments. As will be discussed below, it is the relatively 
low energy per unit angular momentum that is key to the 
proposed metastability of vortex line solutions. 



II. CURVED VORTEX LINE SOLUTIONS 

We now turn to more general vortex line solutions. The 
flow field around a general vortex line in bulk helium is 
homologous to that of the magnetic field around a curved 
wire that follows the vortex line. The "current" in this 
wire is proportional to the quantum of circulation, k — 
h/m. Thus the flow field v(r) is given by the Biot-Savart 
equation 2 ^ 



v(r) 



K 
47T 



(s — r) x ds 



(2) 



The vortex must either form a closed loop or end at a 
boundary of the supcrfluid helium. v(r) must not have a 
normal component at any boundary of the supcrfluid. 
This implies that the vortex must intersect a helium 
boundary at normal incidence. 

A curved vortex will move in its own flow field. When 
the local radius of curvature of the vortex IZ(s) is large 
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TABLE I: Properties of Helium droplets and their straight vortex solutions as a function of the number of Helium atoms. 



He Number 


10^ 


10 s 


10 4 


10 b 


10 b 




Radius R 


1.033 


2.222 


4.785 


10.31 


22.21 


nm 


Helium Binding Energy 


4.72 


6.03 


6.63 


6.91 


7.04 


K 


Total Thermal Ripplon Energy 


0.413 


3.56 


18.5 


87.0 


403. 


K 


L = 2 Ripplon Energy 


1.05 


0.332 


0.105 


0.0332 


0.0105 


K 


Thermal Ripplon y/{L(L + 1)) 


1.54 


8.54 


38.7 


173. 


781. 


h 


Vortex Energy E DFT 


44.3 


129. 


349. 


905. 


2283. 


K 


Vortex Energy E U c (Hollow Core) 


43.4 


122. 


323. 


827. 


2064. 


K 


Vortex Angular Momentum L v 


10 2 


10 3 


10 4 


10 s 


10 6 


h 


E H c/ 'Lv 


56.8 


15.9 


4.23 


1.08 


0.270 


GHz 


Lowest Ripplon E/L 


68.6 


21.8 


6.88 


2.18 


0.688 


GHz 


Vortex-induced eccentricity 


0.53 


0.44 


0.35 


0.27 


0.20 




B/A 


0.85 


0.90 


0.94 


0.96 


0.98 




Deformation Stabilization Energy 


2.00 


3.74 


6.16 


9.22 


12.9 


K 


Max. L loss by 1 atom evaporation 


26 


96 


345 


1198 


4090 


h 


Max. L loss by n atom evaporation 


41 


222 


1213 


6561 


35000 


n 


for evaporation of n = 


5 


10 


25 


59 


146 


atoms 


Max. L loss by fission 


28 


210 


1488 


10150 


67400 


h 


xq for stability limit of curved vortex 


0.72 


0.76 


0.81 


0.85 


0.89 


R(N) 


Min. L for stable vortex 


18 


131 


850 


5.14 x 10 3 


2.96 x 10 4 


h 


Min. E for stable vortex 


9.42 


27.2 


62.1 


126. 


242. 


K 


Max. v for stable vortex 


56.3 


39.0 


27.6 


19.5 


13.8 


m s _1 


Max. kinetic energy of 100 u dopant 


19.1 


9.1 


4.6 


2.3 


1.1 


K 



compared to its core radius a, the motion of the vortex 
core can be calculated using the Local Induction Approx- 
imation (see Appendix iBl and Ref . Ha) . The shape of the 
vortex core surface is determined by the condition of no 
helium flow across it, and turns out to be circular in cross 
section as long as a <C TZ. For finite core size, match- 
ing this boundary condition requires either that the core 
shape be altered or that one add an irrotational solution 
to Laplace's equation that corrects the normal compo- 
nent of velocity. It is not evident to the authors which 
change to make, and so these errors are neglected in the 
rest of the paper. 

We now specialize to the case of a vortex line inside 
a spherical droplet of radius R. Muirhead, Vinen, and 
Donnelly 23 showed that for any arbitrary vortex shape in 
a spherical droplet, the boundary conditions of the flow 
velocity on the surface of the droplet can be satisfied by 
continuing an image vortex outside the droplet. Each 
point on the vortex s (s is the magnitude of this vector) 
generates an image point Si = (R/s) 2 s with vorticity 
equal to —ks/R. Vorticity is conserved by attaching to 
each point along the image vortex a radially pointing 
vortex going to infinity that has a circulation strength 
given by the decrease in circulation of the image along 
its length. These all combine in the Biot-Savart equation 
to give 

K f f (s — r) x ds s (sj — r) x dsi 
V{V) = 4^ i vortox V \s-r\ 3 ' R ' k-r| 3 

s ■ ds r x s \ 

R ' s 2 |si-r| 2 + s|si-r|(i2 2 -r-s)J ' [ > 



where 

dsi = Jj2 ^_2( i _d f ).y (4) 

The three terms in the integral arise from the vortex, im- 
age vortex, and vorticity conserving radial vortex lines, 
respectively, and make contributions that decrease in 
magnitude the order given. 

We seek curved vortex line solutions {x(£), z(£)} that 
rotate at constant angular velocity f2 around the z axis, 
which implies that these vortex line solutions will have 
constant shape. Appendix IBl gives the numerical proce- 
dure used to determine these solutions. Figure [21 shows 
some of the solutions for several values of xo, the distance 
of minimum approach to the z axis. 

Bauer, Donnelly, and Vinen 2 ^ showed that the total 
angular momentum L v and kinetic energy E v of the he- 
lium flow (which are defined by volume integrals for the 
corresponding densities) can be reduced to two surface 
integrals 

L y = mp^n J r x dS (5) 

E v = ~mp b K J v ■ dS, (6) 

where the integration is over the the region in the xz 
plane bounded by the vortex and the surface of the 
droplet, v in Eq. 10 is given by Eq. The origin 

for vector r in Eq. JSJ must be taken as the center of 
the sphere so that the outer surface of the sphere does 
not contribute to this integral expression. These expres- 
sions have neglected a contribution of the integral over 
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the surface of the vortex core, but that should be small 
as long as 1Z ^> a as required by our approximations. In 
the Appendix, we give explicit expressions for the lowest 
order (in a /TV) core surface corrections to L v and E v , 
which are the calculations reported below. For the vor- 
tices considered here, L v is parallel to the z axis owing to 
the symmetry of the vortex and its flow field with respect 
to reflection in the xy plane. Calculation of E v requires 
evaluation of a triple integral, with an integrand that is 
nearly singular along one of the edges of the integration 
domain; one of the integrations can be done analytically, 
leaving a double numerical integration. In the same way, 
the double integral for L v can be reduced to a single nu- 
merical integration. 




x/R 

FIG. 2: Curved vortex solutions for different values of Xo, 
the distance of closest approach to the z axis, in a droplet of 
N = 10000 He atoms. 

The treatment as yet is for vortex lines as classical 
objects with angular momentum pointing in a definite 
direction in space. By the standard rules of semiclassi- 
cal quantization, vortex eigenstates can be constructed 
as linear combinations of the vortex lines with angular 
momentum pointing in all possible directions, with an 
amplitude for each direction given by a spherical har- 
monic. Such states are eigenstates of both the total he- 
lium angular momentum (quantum number L v ) and also 
its projection on the laboratory Z axis (quantum num- 
ber M v ). The vortex lines we have discussed have total 
squared angular momentum that span the range (Nft,) 2 to 
as Xo goes from to R. Semiclassical quantization will 
restrict xq to TV — 1 values with total angular momentum 
quantum number equal to L v = N — 1, N — 2, . . . , 1. The 
straight line vortex is not allowed since the length of the 
total angular momentum of the allowed vortex solutions 
is \L V \ = hyj L V (L V + 1), and the straight vortex solu- 
tion has \L y \ — hN. Each vortex solution with integer 
Ly quantum number has a (2L V + l)-fold M-degeneracy. 



This set of solutions gives N 2 semiclassical vortex states 
for a droplet of N helium atoms. In principle, there are 
other states that involve vibrational excitation of the vor- 
tex lines^i around these Bauer-Donnelly- Vinen vortex 
solutions, but these will not be considered in this work. 

Figure 13 shows the calculated vortex energy and angu- 
lar momentum for a droplet with N = 10 helium atoms. 
It is seen from this figure that the energy of the vortex 
drops as the vortex is moved off axis, going monotonically 
to zero as the vortex is "pushed out" of the droplet. By 
energetic considerations alone, this would imply that the 
vortex solutions are unstable. However, for an isolated 
droplet, one must also conserve angular momentum. A 
vortex can lower its energy by producing a ripplon if the 
derivative of the vortex energy with respect to total an- 
gular momentum, Q' (with units of angular velocity), is 
greater than the E/L ratio of the L = 2 ripplon. Figure^] 
shows Q' for vortices as a function of their minimum ap- 
proach distance to the z axis, normalized to the energy 
per unit angular momentum of this ripplon mode. f2' was 
evaluated by calculation of both E and L for 10 3 values 
of xo for each droplet size, and using Q' — §§-/j^-, with 
each derivative evaluated by finite difference of the cal- 
culated points. It is evident that the vortex is stable to 
ripplon production (has a normalized value of f2' less than 
unity) for most of its range. Tabledlist the minimum en- 
ergy and angular momentum of the curved vortex states 
that are thus stable. 




FIG. 3: Energy (solid line) and angular momentum (dashed 
line) of a curved vortex solution as a function of the distance 
of closest approach of the vortex from the z axis. These re- 
sults are for a, N — 10 4 helium atom droplet. The energy Eo 
of the straight vortex is given in Expression (0; its angular 
momentum is Lo = Nh. 

Because the energy and angular momentum in the hol- 
low core model arise entirely from helium motion, we had 
anticipated that f2' would be equal to 17, the angular ve- 
locity of the curved vortex solution around the z axis. For 
vortex rings in the local induction approximation, Ray- 
field and Reif 18 found that the vortex velocity satisfied 
v = dE/dp, where p is the net helium linear momentum 
due to helium flow around the vortex ring. This insures 
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FIG. 4: Angular velocity as a function of distance of closest 
approach, calculated for droplets with N = 10 2 , 10 3 , 10 4 , 10 5 , 
10 6 (bottom to top) He atoms. Plotted is ujR 2 = 4tvQR 2 /k, a 
dimensionless quantity. Solid lines are computations using the 
model described in Appendix [B] dotted lines use the model 
of Ref. I22I As discussed in the Appendix, the local-induction 
approximation breaks down as xq —* R. 
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FIG. 5: Comparison of the vortex "angular velocity" fl' — 
AE/AL and the angular velocities ujl of the lowest ripplons 
(L = 2 . . . 7), as a function of distance of closest approach to 
the 2 axis, for droplets with N = 10 2 , 10 3 , 10 4 , 10 5 , 10 6 , 10 7 
(top to bottom) He atoms. All angular velocities are shown 
as ratios to lo 2 = \/ 8ira / (3mN) . 



that as xo — > R (but R — xo ^> a so the local induction 
approximation still holds), Q' — > Q. For a vertical vor- 
tex in a cylinder (where E and L have simple analytical 
solutions), we have demonstrated that Q' = Cl'. For our 
solutions of curved vortex lines in a spherical droplet, we 
however find that O' ^ fl', as demonstrated in Figure 
where both are compared as a function of xo / R for var- 
ious droplet sizes. It appears from our calculations that 
Q! — > £1 as N — > oo (except for the log divergence of il as 
xo — * 0), but the rate of convergence is rather slow. It re- 
mains to be established whether the discrepancy between 
f2 and f2' is due to an error in calculation or whether it 
is real. 
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FIG. 6: Comparison of the computed reduced angular veloc- 
ity Lj — i^fi (dashed lines) and the quantity u>' = ^-AE/AL 
(solid lines), for droplets with N = 10 2 , 10 4 , 10 6 atoms (bot- 
tom to top). 

It is evident from figure that f2 peaks for the straight 
vortex (xq — 0); more careful analysis shows that there 
is in fact a logarithmic divergence: ft sa „ K pa In as 
xq — * 0. In the limit that the vortex line moves near the 
surface of the droplet, our solutions no longer provide 
meaningful estimates as the vortex ring radius R — xq 
approaches the value of the core radius a. In addition to 
the expected breakdown of the local induction approx- 
imation in this limit, such a vortex would be localized 
in the region of highly inhomogeneous density near the 
helium surface. Feynman has suggested that the limit of 
a vortex ring of atomic scale is the roton quasiparticle 
excitationiS, 

In order to test that there was no numerical error in 
the evaluation of the integral expressions for L v and E v 
reported in this work, calculations for selected values of 
a/R were performed with totally independent routines 
written by the two authors, using distinct programing 
packages (Mathematica and Mathcad), and the values 
of fi, L v , and E v were found to agree to six significant 
figures. 

III. POSSIBLE SPECTROSCOPIC 
SIGNATURES OF MOLECULES BOUND TO 
VORTICES 

It is well know in bulk liquid helium that impurities 
tend to bind to vortex lines or vortex rings. The bind- 
ing of ions to vortex rings was the basis of the remark- 
able experiment of Rayfield and Reifi^ that measured 
many of the properties of ring vortices. Calculations by 
both Dalfovo et alA and by Draeger and Ceperlej*^ have 
demonstrated that atomic or linear molecular impurities 
in a helium nanodroplet bind to linear, maximum an- 
gular momentum vortex lines. For technical reasons, 
both these calculations required cylindrical symmetry, 
and thus could not study the energy change upon ro- 
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tation of the molecule with respect to the vortex axis. 
However, it is generally believed that the binding energy 
is maximized by aligning the linear molecular axis with 
the vortex since this displaces the maximum amount of 
helium from near the vortex core, where it has the high- 
est kinetic energy; however, the increased helium density 
in the first solvation layer around the solute may in part 
counteract this effect. This suggests that the anisotropy 
in the binding energy should be of comparable magni- 
tude as the binding itself, but this natural expectation 
has yet to be checked by calculations with a realistic he- 
lium density profile. It seems likely that the magnitude 
of the anisotropy of the binding for linear molecules, such 
as HCN and HCCCN, is much higher than the rotational 
constant of these molecules in liquid helium. If this was 
correct, then the gas phase- like rotational structure of 
these molecules would be quenched, leading to a rovi- 
brational spectrum dominated by a Q branch; such Q 
branches are, however, absent in the observed spectra 
of these molecules i^i This suggests that at most a small 
fraction of the droplets probed in these experiments could 
have such a linear vortex line excitation. It is further 
noted that the vibrational spectrum of the spherical top 
molecule SF 6 is also expected to be a sensitive probe for 
the presence of a vortex line, as the highly anisotropic 
helium density created by the vortex core would lift the 
triple degeneracy of IR active fundamental modes, pro- 
ducing a spectroscopic structure in qualitative disagree- 
ment with what has been observed^ 

Molecules are also expected to bind and align with 
curved vortex lines, as long as the solute kinetic energy 
required for it to move with the curved vortex line does 
not exceed the binding energy of the solute molecule to 
the stationary vortex. The velocity of a curved vortex 
increases smoothly from zero when xq = to a maxi- 
mum value (given in table [1} when the vortex becomes 
unstable to ripplon formation. For a molecule with a 
translational mass of 100 u, the maximum kinetic energy 
that the molecule requires to stay bound to a curved vor- 
tex line is between 19.1 and 1.1 K as N varies between 
10 2 and 10 6 (see Table HJ. This kinetic energy can be 
compared to binding energies to the straight vortex of 
5.0, 4.4, and 7.7 K for Xe, HCN, and SF 6 , respectively, 
as calculated by DFT. 4 Thus, for the least stable curved 
vortex lines in smaller droplets, the molecules may be- 
come unpinned. This suggests that it would be useful to 
have improved binding energy estimates for molecules, 
particularly to curved vortex states. 



IV. VORTEX FORMATION MECHANISMS 

A possible explanation for the failure of existing exper- 
iments to detect droplets with vortex lines is that they 
cannot be formed in the first place, perhaps because of 
the high energy and angular momentum required. We 
will consider the probability of events that deposit suffi- 
cient energy and angular momentum in the droplets such 



that vortex formation can in principle take place. 

In most experiments, the droplets grow out of a super- 
cooled gas^&In this situation, the large droplets probably 
grow at least in part by coalescence of smaller droplets. 
If we consider the coalescence of two droplets with N/2 
atoms that collide with an average impact parameter 
(2/3 of either's diameter), the minimum relative veloc- 
ity required to deposit an angular momentum of NTi in 
the combined droplet is v Te \ — 270-/V~ 1 / 3 m/s, which is 
58 - 2.7m/s for N = 10 2 — 10 6 . Experiments of Toennies 
et almL have found that the expansions produce droplets 
with a final speed ratio of rs 100, which implies final 
relative velocities of w 8 m/s for expansion from a 20 K 
source. While such relative collision velocities are, except 
for the largest droplets, less than those required to pro- 
duce a vortex line with maximum angular momentum 
(straight vortex), they are sufficient to produce curved 
vortices. Further, the above estimates are only average 
values, and the relative velocities in the part of the ex- 
pansion where the droplets undergo substantial growth 
are likely significantly higher. Experiments that produce 
droplets above N w 5 x 10 4 typically use expansion con- 
ditions such that liquid helium is ejected through a cold 
nozzle into vacuum, and fragments by cavitation as the 
pressure falls far below the equilibrium vapor pressure. 
This break-up is believed to generate considerable turbu- 
lence, and thus may lead to dense vortex formationi^SiSa 

The "pick-up" process^ by which droplets are doped 
with solutes, should often deposit enough angular mo- 
mentum, and almost always enough energy, to form a 
linear vortex line, particularly for the pickup of heavy 
molecules by not too large droplets. It has been pre- 
dicted that vortices are nucleated when the velocity of 
an impurity exceeds to sound velocity in helium, 31 which 
is well below at least the typical impact speed of atoms 
or molecules striking helium nanodroplets. Consider the 
pickup of atoms or molecules of mass M from a thermal 
gas at temperature T by droplets with laboratory frame 
velocity v&. Integration over the Maxwell-Boltzmann 
distribution for the gas gives the probability density P r 
for the relative velocity v r between droplet and gas 



Pr(Vr) 



2M v T 
nkT Vd 



sinh 



MVrVd 




kT 


exp 



2kT 



_ . (7) 

From this distribution, the average relative collisional 
speed, v T , is calculated to be 



2kT 



exp 



Mvl 



2kT 



+ 



kT \ , Mv 2 A 

Vd+ Mv- d ) eii y^kT- 



For a given impact parameter b, all collisions with rel- 
ative velocity v r > v m (b) = %/b ■ (1/m + N/M) result 
in angular momentum of at least NTi, enough to form a 
straight vortex. The fraction of such collisions is 



*(6) = 



JZ(b) VrPr(v r )dv r 
/ °° V T Pr(v T )dv T 



(9) 
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If we assume that the the pickup probability is inde- 
pendent of collisional impact parameter b for b < R = 
tqN 1 / 3 , and falls abruptly to zero for b > R, then the 
fraction of resulting doped droplets that have at least 
Nh of angular momentum from the pickup process is 



>Nh 



R 2irbdb 



$(6). 



(10) 



Figure [7| shows plots of the fraction of pickup collisions 
that result in an angular momentum greater than Nh, as 
a function of N for a number of solutes. Figure |S| shows, 
for droplets of N = 10 4 helium atoms, the distribution 
of collisional angular momenta deposited by pickup of 
the same solutes. In many cases relevant to previously 
reported experiments, the pickup process has a fair prob- 
ability to deposit sufficient angular momentum to create 
a straight vortex, and almost always sufficient angular 
momentum to create curved vortex lines. Thus, the lack 
of observation of vortex lines does not appear to be due 
to a lack of initial angular momentum to form them. 




1000 

number N of He atoms in droplet 



10000 



FIG. 7: Probability of a pickup collision depositing more than 
Nh units of angular momentum in the droplet, as function of 
droplet size N. Droplet velocity = 456 m/s (20 K super- 
sonic expansion); gas temperature T — 300 K. 



V. STABILITY ANALYSIS OF VORTEX STATES 

An alternative explanation for the lack of evidence of 
vortex lines in experiments to date could be that vortex 
lines are unstable, and decay or are expelled from the 
droplets in a time short compared to the time between 
the pickup of solutes and the spectroscopic detection of 
the doped droplets. In this section, several possible vor- 
tex destruction pathways are considered, concentrating 
on the simplest case of the linear vortex line. 

A helium nanodroplet excited with a linear vortex line 
has an angular momentum of Nh, much higher than the 
magnitude of the angular momentum thermally present 
in the ripplons at the droplet temperature of ~ 0.38 K 
that has been experimentally found for helium nan- 
odroplets cooled by evaporation. 6 For droplets in the size 
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FIG. 8: Probability densities of collisional angular momenta 
during pickup of various molecules by a droplet with N = 10 4 
He atoms. Droplet velocity v$ = 456 m/s (20 K supersonic 
expansion); gas temperature T = 300 K. For other droplet 
sizes, the shapes of the distributions are the same except for 
a correction due to changes in the reduced mass, and scaling 
of the average collisional angular momentum as N 1 ^ 3 . 



range of experiments reported to date, the thermal exci- 
tation of phonon states is negligible in comparison to that 
of the ripplons^S The lowest frequency (L = 2) ripplon 
mode has the lowest energy per unit angular momentum 
of any of the helium excitations except vortices, but this 
value is still higher than the corresponding value for the 
linear vortex, as indicated in Figure □ Tabic [Q shows the 
ratio of excitation energies and angular momenta for the 
L = 2 ripplon and the linear vortex for a range of helium 
droplet sizes. While decay of the vortex state is highly 
exothermic, this is only true if one ignores the require- 
ment of conservation of angular momentum; this neglect 
is appropriate for helium contained in a vessel with walls, 
to which the angular momentum can be transferred, but 
not for a droplet isolated in a high vacuum chamber. 
When angular momentum conservation is enforced, the 
decay of the linear vortex line into ripplons is highly en- 
dothermic and thus cannot occur for a droplet unless 
the droplet also contains a counter-angular momentum 
in ripplons that is nearly as large as the angular momen- 
tum of the vortex itself. Table[I]also contains the thermal 
average value y/ (L{L + 1)) of ripplon angular momentum 
for different size droplets at 0.38 it is only a small 
fraction of that of the straight vortex. Exchange of an- 
gular momentum between thermally populated ripplons 
and a vortex could at most allow the vortex to move 
slightly off axis. 

Helium nanodroplets can shed angular momentum by 
helium atom evaporation. A helium atom evaporating 
with momentum p can carry away a maximum angular 
momentum of rs pR, which reduces the internal energy 
of the droplet by E\, + p 2 /(2mn e ), where E\, is the bind- 
ing energy of the helium atoms to the droplet. For finite 
droplets, Ey, is reduced from the bulk value due to a term 
that reflects the change in surface energy upon evapora- 
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tion of an atom. 37 Table [I] shows the maximum angular 
momentum that can be lost by single atom evaporation 
using the entire energy of the linear vortex line^S In every 
case, that angular momentum is far lower than the angu- 
lar momentum of the vortex line that must be carried off 
for the vortex line to decay. Evaporation of multiple he- 
lium atoms will increase the angular momentum that can 
be carried off for the same total helium kinetic energy. 
Because the available kinetic energy is reduced by the 
increased binding energy of the multiple helium atoms, 
there is a number of evaporated atoms that allows the 
maximum angular momentum to be removed. Table [I] 
also lists this amount of angular momentum for a range 
of droplet sizes M> It is evident that in all cases, this max- 
imal amount of angular momentum lost by evaporation 
is far below that which is needed for the vortex line to 
decay while conserving angular momentum. 

The above analysis assumed that the helium atoms 
evaporated as isolated atoms. An alternative decay 
mechanism is fission of the droplet: in this case, the 
loss of helium binding energy is only due to the increase 
in surface energy, which is higher for the two droplets 
than for the original nearly-spherical droplet. Fission of 
a droplet into two equal-size droplets costs far more en- 
ergy than that contained in a vortex line, even neglecting 
the required kinetic energy for relative motion of the frag- 
ments. Table [I] contains the fragmentation that carries 
away the maximum possible quantity of angular momen- 
tum, taking the original droplet radius as the impact pa- 
rameter of the departing fragments^ Yet again, this is 
far below the angular momentum of the linear vortex. 

The above considerations show that the simultane- 
ous constraints of conservation of energy and angular 
momentum prevent decay of straight vortices in a pure 
droplet by any of the mechanisms that are known to the 
authors. Adding a molecule to the droplet will further 
decrease the energy per unit angular momentum of the 
vortex as long as the molecule remains bound to the vor- 
tex. It thus appears that droplets with a vortex line 



should be stable to decay. 

How can we understand the failure to date to ob- 
serve vortex lines in helium nanodroplet isolation spec- 
troscopy? The answer is probably related to the fact that 
under typical experimental conditions, the above mecha- 
nisms, by putting sufficient angular momentum into the 
droplets to form vortex lines, add far more than the min- 
imum energy required to do so. The density of ripplon 
states grows very rapidly with increasing energyjS it may 
well be that at the energy and angular momentum of a 
droplet following coalescence or pickup, the fraction of 
states that contain a vortex line is a totally insignificant 
fraction of the total density of states, and that droplets 
will shed sufficient angular momentum along with energy 
as they evaporatively cool that they have a negligible 
chance to be trapped in one of the metastablc solutions 
containing a vortex. We have recently carried out Monte 
Carlo cooling calculations, using statistical reaction rate 
theory, for the evaporative cooling of pure and doped he- 
lium nanodroplets, conserving both energy and angular 
momentum's considering only ripplon excitations of the 
droplets. In the future, we hope to extend that work to 
include the spectrum of vortex line excitations as well, 
which it is hoped will shed further light on the still un- 
solved problem of why spectra of molecules bound to vor- 
tex lines has yet to be observed in Helium Nanodroplet 
Isolation Spectroscopy. 



APPENDIX A: HOLLOW-CORE MODEL 

In the local-induction approximation, the flow field 
around a curved vortex is approximated locally by the 
flow around a vortex ring of equal radius of curvature. 
The velocity potential for the flow around a vortex ring 
in the x — y plane, with radius Tl, centered at the origin 
and moving in the +z direction, can be computed from 
Eq. ©: 



nTlz 



tt(K 2 -r 2 - z 2 )y/{Tl + r) 2 +2 S 



+ (n - vV + z 2 ) ■ s 



+ z 2 )-E 
(Tl 



(Tl - rf + z 2 



Tl 2 + r 2 + z 2 
- r) 2 + z 2 



2-rVT 2 



(Tl - r) 2 + z 1 
1*' (Tl + r) 2 + z 2 



(Tl - r) 2 + z 2 
K 2 +r 2 + z 2 + 2'llV7 r +1?' {K + r) 2 + z 2 



-sign(z), (Al) 



with the function 



(r,z) = Tl ■ (1 + Ccostf,Csintf) and < C < 1, yields 



n(n|m) 



' n(H|I) 

m mm 



(A2) 



in terms of complete elliptic integrals of the third kind, 
n(n|m); k = h/m is the quantum of circulation. Series 
expansion of Eq. (|A1|) around the vortex singularity, with 



«0 <, 8 

In — sin if 

2tt 4tt C 



^(ln^-§)sin2tf. (A3) 



In the hollow-core model, we assume that the vortex sin- 
gularity is surrounded by an empty "vortex core" of ra- 
dius a, outside of which the helium density is constant 
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and equal to the bulk value, pt>- The helium velocity 
field is given by v = —~V(j>: the first term in Eq. I|A3[) is 
the source of the looping helium flow around the vortex, 
the second term gives rise to the forward motion of the 
vortex ring with speed 



a 



(A4) 



along the +z axis, and the third term allows us to es- 
timate the residual flow across the vortex core surface 
due to the fact that the exact vortex core surface is not 
circular in cross section: the ratio of the rms velocity 
of this residual normal flow to the forward velocity is 
3a/{4TlV2) as a/K -> 0. 

Rayfield and Rei£i& have experimentally determined 
the forward velocity u of ring vortices, and determined 
a core radius of a = 1.00 A by using the formula 
u = [hi — |] . Comparison of this expression to 
Eq. (|A4|) yields a = a/^fe = 0.607 A, which is the core 
radius that we use in this work for the hollow-core model. 



APPENDIX B: VORTEX SHAPE CALCULATION 



pick the parameter u>, and thus the angular velocity $1, 
such that the resulting vortex line intersects the droplet 
surface perpendicularly; for a spherical droplet with ra- 
dius R, this condition is x{£i)z'{£i) = z(£i)x'(£i) with 
x(£ 1 ) 2 + z(£ 1 ) 2 = R 2 . 

This procedure differs from the one given in Ref. |52 
in that Bauer et al. neglect the variation of the Lambert 
function, and replace the denominator of Eq. i|B3(l by a 
constant. For droplets consisting of 100 < N < 100000 
helium atoms, the error of this approximation is below 
0.6% for the evaluation of the angular velocity. 

Figure 0] shows angular velocities of curved vortices 
computed with both models. 

As xo — * R, the vortex solutions are close to half-circles 
with radius R — xq. We expect the hollow-core model to 
break down if R — xq < a. A clear sign of this breakdown 
is that for xo — > R, the Lambert function in Eq. I|B3|) 
starts giving complex values as the droplet surface is ap- 
proached (aujx(£)e/8 > 1/e). In the simplified model of 
Ref. the angular velocity turns toward negative val- 
ues for R — xq < a, invalidating that model as well. It is 
to be noted that in this limit, the description of the vor- 
tex core as a cylindrical tube around the vortex is very 
inaccurate. 



A vortex, s(£), is most generally described as a curve 
in the x, z plane parametrized by its length £; three- 
dimensional vortex curves are always longer and thus 
higher in energy for given boundary conditions^* The 
local vortex curvature is 



xW 



x"{£) _ z"{£) 
z'{£) ~ x'(£)' 



(Bl) 



In the hollow-core model with the local-induction approx- 
imation, the magnitude of the vortex velocity is deter- 
mined from the core radius and the local curvature, and 
its direction is given by the binormal vector (the normal- 
ized cross product of ^| and ^f). The magnitude of the 
velocity is given by Eq. JHJ|: 



In 
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(B2) 



For a vortex to rotate around the z axis at angular ve- 
locity f2 without changing shape, it must satisfy v{£) = 
flx(£), which leads to the coupled differential equations 
that determine the vortex shape: 



{x"(£),z"(£)} 



ux{£) 



r{-z'(£),x'(£)}, 



-awx(t)e/%] 

(B3) 

where we have used uj = A-kVI/k, and the Lambert func- 
tion W-i(x)^ defined as the smaller of the two real 
solutions of x = ye y for —1/e < x < 0. Combined 
with the initial conditions {x(0), z(0)} = {xo,0} and 
{.t'(0), z'(0)} = {0, 1}, Eq. CE| yields the desired vortex 
shape functions by numerical integration. For a given 
minimum approach distance xq to the z axis, we must 



APPENDIX C: CORE CORRECTIONS 

Bauer et al. show— that the angular momentum and 
kinetic energy of a droplet with vortex can be written as 
surface integrals: 



L v = mp\> I (f>r x dS 
Js 

E v = ^mp h J 4>v ■ dS, 



(CI) 
(C2) 



where E is the surface of a connected volume of helium 
with no branch cuts in the velocity potential 4>, and v = 
— V0. While most of the angular momentum and kinetic 
energy comes from the surface integrations of Eqs. (|5I6|I . 
there are corrections due to the surface integrations of 
Eqs. ljCHC2jl over the vortex core. Using Eqs. (|A^|) . (JbTJ, 
and 1B3I). we find these corrections to be 





mpK 




auj 2 


SE V p 


mpK 2 




8irauj 



Ci,o 



^2,0 + -jp 21 



with 

C„, m = {auj) n+l 



"Cl,0 + ^^2,0 — g^ 2 ' 1 



x{£) n 
WH\(-aLux(£)e/8) 



(C3) 
, (C4) 

d£. (C5) 



From symmetry considerations, we must have 
L v , z (— x ) = L VjZ (x ) and E v (-x ) = E v (x a ). In fact, 
the first core correction (C^q) adjusts L v z to satisfy 
this symmetry constraint, i.e., L v , z (xq) 



Nh - 0{x 2 ) 



10 



However, this symmetry is not exactly satisfied for work was supported by a grant from the National Sci- 
E v (xq), which we assume is due to inaccuracies of the ence Foundation, 
local-induction approximation. 
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